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Abstract 

Applying a bi-local mean field approximation to the fractional quantum Hall state of 
V = 1/3, we obtain charged and neutral vortex mean field solutions numerically. We calculate 
the mean field energy and the fluctuation corrections. The charged vortex has a fractional 
charge and a fractional angular momentum. The neutral vortex is a bound state of two 
charged vortices and has a zero charge and a zero angular momentum. The creation energy 
of the neutral vortex is about a half of the pair creation energy of two charged vortices. The 
magnetic field dependence of the gap energy agrees with the Laughhn's quasiparticle gap 
energy. 
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1. Introduction 



We consider the two-dimensional interacting spinless electron system in a perpendicular 
magnetic field. The fractional quantum Hall effect (FQHE) can be characterized by the 
existence of an energy gap in a partially filled Landau level. Strong correlations among 
electrons with the Coulomb interaction are responsible for the energy gap. Since a par- 
tially filled many-electron state in the Landau level has an enormous degeneracy, systematic 
calculations based on a microscopic theory are very difficult. Laughlin proposed a phe- 
nomenological many-body wave function for a ground state of an incompressible quantum 
liquid^\ In the Laughlin's theory the excited state includes quasiholes and quasiparticles 
which have fractional charges and obey the fractional statistics. The pair creation energy 
of a quasihole and a quasiparticle gives a finite energy gap and is in a good agreement 
with experimental values at strong magnetic fields. Although the Laughlin's theory explains 
the FQHE phenomenologically, it may be desirable to have more systematic calculational 
method. 

Recently a bi-local mean field theory was applied to the many-body system for the 
FQHE and obtained a uniform mean field solution for a ground state'^'' and charged vortex 
solutions for a topological excited state^^. It was found that the vortex solutions have a 
fractional charge and a fractional angular momentum and have vorticity one. In this theory 
the correlation function {il)"^ [y)il)[x)) is treated as a mean field. It was found that the gap 
energy is rather smaller in a mean field theory than the Laughlin's one. In the present paper 
we calculate the fiuctuation corrections to the vortex energy and find that the gap energy 
becomes close to the Laughlin's one. We also calculate charged and neutral vortex solutions 
of vorticity 2. The neutral vortex solution has a zero charge and zero angular momentum. 
The neutral vortex is considered to be a bound state of two charged vortices which have 
charges of opposite sign. The neutral vortex energy is about a half of the pair creation 
energy of two charged vortices. 
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In §2, we review a bi-local mean field theory for the FQHE. We point out that we need 
to change the self-consistent condition for a systematic expansion around the vortex mean 
field solutions slightly. Nevertheless, it is shown in appendix that our results obtained in 
this paper is not changed. In §3, we consider the ansatz forms of charged vortex solutions of 
vorticity riy and get the exact relation between the fractional charge and fractional angular 
momentum of vortex in the large B limit. We also consider a neutral vortex of vorticity 
2. The fluctuation corrections to the vortex energy are calculated in §4. The results of 
numerical calculations are presented in §5. In §6, we summarize our results. 

§2. Mean field approximation for ttie FQHE 

We consider the two-dimensional electron system in a perpendicular uniform magnetic 
field B which is described by a Hamiltonian 

V) = J <fx[^Hx)^-^^^^i;{x) + (fyi;\x)i;\y)V{x - (2.1) 

where Pi — A — ^{—X2, xi) and V{r) — e^/ («;r). We use the unit h — c — 1 and fix 

K = 13, m = 0.07me for GaAs. We supposed that the magnetic field B is strong enough to 
make electrons fully polarized and neglected the spin effect of electrons. 

In the functional integral formalism, the partition function is written as 

Z = Tr(e-^^(^^'^)) = J Vi^^Vi^Q-K^[- dt{ j d^xi^^dti^ + H{i^\^))], (2.2) 

where tf^ and if) are anti-commuting field variables, ■0^ and are anti-commuting field 

operators and j3 = l/ksT. We introduce a bi-local auxiliary field U{x,y]t) and convert the 
interaction term into a quadratic form of -0^ and -0 as 

Hu{i^\i^) = j <fx^\x-t)^^^^^^^{x-t) 

+ 11 £x(fyV{x - y) [f/(x, y- t)U{y, x; t) ^^-^^ 
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- U{x, y; t)^^{x; t)^{y; t) - U{y, x; t)^p\y- t)^{x; t)], 



Z = Af-^ J ViP^VipVUexpl- dt{ J (fx^p^dtip + Hu{ip\ ^))] 

j V^^'D^VUex.p[-S{ip\^,U)]^N'-^ J VUexpl-SeffiU)], (2.4) 

Af^ J VUex^[- j\t{ j d^xd^y ^^"^ ~ U{x, y; t)U{y, x; t))] . (2.5) 

In the mean field approximation, the partition function is calculated at the stationary point 
of S^ff as 

Zo = e-^««(^°) = Tr{exp[-pHuS^,m, (2-6) 
dS,s 



= 0. (2.7) 

U^Uo 



dU 

Using a correlation function, Eq.(2.7) can be written as 

Uo{x,y) = {ilj\y)ij{x))u„ (2.8) 

^Hymx))^^ = Zq' J V^l;^V^l;{^l;\y)^{x))exp[-SiiJ\iJ,Uo)] 

= Z^'Tr{^\y)^ix)exp[-pHuS\m {Eo\tl'HyMx)\Eo) , (2.9) 

in the hmit /3 — > oo. The state \Eq) is the lowest energy eigenstate of Huo{i^\ ■0). We solve 
the self-consistency condition Eq.(2.8) under the following ansatz form : 

C/o(x, y) = Uqp{x, y)e-''^^-y^" exp[i ^^(0^^], (2.10) 

Jx 

where Uq — u/ttRq, 7 = 1/2-Rq, Rq — ^/2/eB and p{x,y) — p{y,x). The line integral is 
along a straight line between x and y. p{x,y) and are unknown real functions and we 
will determine these functions under appropriate assumptions. We assume that the mean 
field UQ{x,y) coincides with a uniform self-consistent solution at infinity : 

lim p(x, x) — 1, 
hm aiiO^eMO- ^ ' ' 

Equations (2.10) and (2.11) mean that the electron density Uq{x^ x) becomes uniform density 
of the filling factor v at infinity. By substituting Eq.(2.10) into Eq.(2.3), we get a mean field 



Hamiltonian : 

Huo^ J (fx^Hx)ho^{x) + ^ J d'xd^yV{x-y)p\x,y)e-''^^--y^\ 
ho = [^^-^^ - + af)p{x, (2.12) 

3/2 

where Iq is a zero-th order modified Bessel function, ho is a single-particle Hamiltonian and 

[ ]x^x represents to take a limit x ^ x after P operates on a function of x. 

It is convenient to expand the electron field operator by the /iq's eigenfunctions. We 
assume that the system has a rotational invariance and the eigenfunctions are classified by 
the angular momentum quantum number, I, as 



l,n 

where dj^ and di'^n' ^-re anti-commuting creation and annihilation operators. The N-electron 
state is defined by a totally anti-symmetric tensor -F;i,.../jv ^ 

M 

E Fh,-iAAfi- --^Loi^)^ ..... 

ai,n |0) = for any I, n, 

where M is a maximum value of the angular momentum of a single-particle. Since we 
concentrate upon the case oi u < 1 (especially u — 1/3), we restrict \^) within the n — 
eigenstates in Eq.(2.14). For a free electron system the n = eigenstates belong to the 
lowest Landau level. .../^ is zero unless h + h + ■ ■ ■ + In — Lt for a fixed total angular 



momentum Lt- Prom this property one can prove the following relations : 

M 

(TV. TV!) J2 Kw.-lMPrnM,---lN - ^nSn^rn, < < 1, 

/2-/^=0 .2.15) 

1=0 

We don't need to know the explicit form of Fi^^...i^ but need to know the value of ui for the 
following calculations. By Eqs.(2.14) and (2.15) we can calculate the correlation function as 

M 

{^\i,\y)i,{x) \^) = J2uiulo{y)ui,o{x). (2.16) 

1=0 

The condition that the electron density becomes a uniform density of filling factor u means 
lim ui ^ u ^ N/(M + 1) for ^ oo. Using Eqs.(2.10) and (2.16), the self-consistency 

l—^O0 

condition (2.8) becomes 

ry ^ 

C/op(x,y)e-^(--2')'exp[z / ai{Ode] = Y.mloiy>i,o{^)- (2-17) 

•^^ 1=0 

Generally the N-electron state |^) in Eq.(2.14) is a superposition of different m+iCn 
states. Hence the different states in l^*) have to be degenerate in the energy eigenvalue be- 
cause the lowest energy eigenstate survives in the self-consistency condition (2.8) as Eq.(2.9). 
We will show later that the energy eigenvalue Ei q is degenerate for a large Z, but is not de- 
generate for a small / in the vortex mean field solution. Therefore we need to modify the 
self-consistency condition (2.8) so as to make the eigenvalues degenerate. In fact this is 
possible by using an arbitrariness in the limiting procedure of discretized time slices^) . The 
details are summarized in appendix. Since the results obtained by using Eq.(2.8) become 
the same results as obtained by using a modified one, we continue to use Eq.(2.8) in the 
following. 
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§3. Vortex mean field solutions 

Following Ref.5), we work on the ansatz form of a vortex as 

ai{x) = eAi{x) + ni{x), n{x) = —^{-X2, xi), (3.1) 
with the boundary conditions : 

n(0) = ny = 1, 2, 3, ■ ■ ■ , 

lim n{r) = 0, (3.2) 

r— >oo 

p{x,0) = 0, p{x,x) = p{r). 
In Eqs.(3.1) and (3.2), the form oini{x) is a reminiscent of the vortex solution in the Abelian 
Higgs theory It is seen by separating the non-singular part fii in rii as 

Jx 

~ , \ "^{^^ ~ ^(0) / s 

^»(^) = \.2 (~^2, xi), (3.3) 



X 



r=oo 



The singular part of rii gives a phase dependence e*"''^ to the expectation value of the field 
operator and the non-singular part gives a dynamical magnetic field, | V x n| = localized 
around a vortex. Uy is called vorticity and characterize the topological excitations. 

The uniform density mean field solution was already obtained in Ref.3) as 

p{x,y) = 1, 

ai{x) = eAi{x), (3.4) 

1^1 — ly for Z = 0, 1, • • • . 
We denote the uniform density state as 

M 

l*o)- E <^..z.«lo<o---4.,olO)' 
The explicit form of F^^^ ^^^ is obtained from the coefficient _ of the Laughlin wave 



function as Y[i<j(^i ~ ^j)^ — Ylih-i^ -In^^^ ■ ■ ■ z^'^ for odd integer P — The eigen- 
function of for the uniform density mean field is written as 



(0) 



^o(x) = (7rit:^Z!)-VV^e-^*^/2, Z = 0,1,2,..-, 

2; = (xi + ix2)/-Ro, 



(3.6) 



and the energy eigenvalue is degenerate as 



^ht-^^^^^- (3.7) 



(A) Charged vortex 

Using the uniform density state we construct a quasihole state of a charge —n^ve as 

M 

\^-n^v)^ XI ^zi'!--«jv^Zi+nv,0^Z2+nv,0 • • • ^L+nv,0 |0) • (3-8) 
1\,Ii---In=0 

Note that in Eq.(3.5) and in Eq.(3.8) are not same because is dependent upon the 
set of eigenf unctions of /io, {'Uz,n(3^)}) ^-^id /iq is dependent upon the state of electrons. The 
state gives vi^ — ■ ■ • — v^^-x — and v^^ — i^n^+i — ■ ■ ■ — y ■ In fact the charge 

difference from the uniform density state is 



nv — 1 

1=0 

In the large B hmit, the single-particle state are restricted in the lowest Landau level and 

;o) 

',0 



Sxei)^ \^-n.v) - (^'ol / Sxei}'^i) l^'o) = e ^ {-v) = -n^ue. (3.9) 



the eigenstate becomes uf'}{x) of Eq.(3.6). Thus the correlation function for a quasihole 



state becomes 



/""^ (3.10) 
= ^(1 - E ^-^e-'y<) exp[-|z, -zy\' + ij eAdC], 



as S — > oo. Prom Eq.(3.10), p^{x,y) reads 



L(x) = tijXiPjiyx). 

Using Eqs.(3.11) and (3.12), the angular momentum of the vortex is given by 



(3.11) 



p\x,y) ^1 - 2e-^^^^-^^ ^^^os{W - Cx^ySine) 

1=0 

+ g-2C.e.cose y ikkL_cos(/-m)^, 
^ l\m\ ^ ' 

l,m=0 

as S — > oo, where Zx — Ck^*^^, 9 — 9x — Oy. For \^x — Cyl ^ reads 

e^"V(nv-l)! 
n(r) ^ — — —, as B ^ oo, 

e^^ - E e^'A! (3.12) 
z=o 

n(0) = ny. 

The momentum and the angular momentum density operator are defined by 



(3.13) 



(^-nvz^l j dPxL{x) \^-nyy) = -Uq J dx'^p{r)n{r) -UyU, (3.14) 
as S — > oo, where we used a relation i^ip^ {x)-^il){x)) = lim ^ {^p\y)ip{x)Y 



Similarly we construct a quasiparticle state of a charge +nve for i/ = 1/3 as 

M 

|*+nvi^) = flriv.O ■ ■ ■ ^2nv-l,0 XI -^Zi,-^]v^/i+2nv,0^/2+2nv,0 ' ' ' ^Zjv+2nv,0 1^) " (3.15) 

h,l'2---lN=Q 

The state \^+n.v) gives z/q = ■ ■ ■ = z/„^_i = and Vn^ = ■ ■ ■ = z/2„,-i = 1 and i'2n. = 
^2nv+i = • • • = = 1/3. In fact the charge difference from the uniform density state is 

/„ riv— 1 2nv— 1 

(fxe'^^ \^-n^^) - (*o| / <fxe'^^ |*o) = e{ ^ (-i/) + ^ (1 - i/)} = +71^1/6. 

Z=0 nv 

(3.16) 

In Eq.(3.16), it is essential that v is equal to 1/3. In the large B limit the correlation function 



for a quasiparticle state becomes 



l—Uv 



(3.17) 



as S — > oo, where 



1, for < Z < riv - 1, 
{ -2, for ny<l<2ny- 1, 
0, otherwise. 



From Eq.(3.17), p^{x,y) reads 



z=o 

Z,m=0 



(3.18) 



as S — > oo. For — Cj/I ^ 1; '^('^) reads 

3e^"V(nv-l)!-2e^"V(2nv-l)! 
n(r) ^^^-^ , as B ^ oo, 

e^'- E c^^e'V^! (3.19) 
z=o 

n(0) = riy. 

Using Eqs.(3.18) and (3.19), the angular momentum of the vortex is given by 



{^+n.u\ J d'xL{x) = -Uo J dx'p{r)n{r) ^ +n^iy, (3.20) 

as B ^ oo. 

(B) Neutral vortex 

We also construct a neutral vortex state for i/ = 1/3 as 

M 

l^neu) = 4,0 E <...Z.«l,0«L3,0 " " " «l+3,0 1°) ■ (3-21) 

The state |^neu) gives uq — ui — 0, U2 — 1, us — U4 — • ■ • — ly — 1/3. The charge difference 
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from the uniform density state is zero. In the large B hmit, p^{x, y) reads 

p\x, y) ^1 - 2e-^^^y ^{008(^x^2/ sin 9) + i^iy cos{e - CxC?/ sin 9) 

- iUyf 008(26 - Uy^nd)} + e-2^^^^^°^^{l + (Uy? (3-22) 
+ (CxCy)^ + "^^x^y COS 9 - 2{^^^yf cos 9 - 2{^^^yf cos 29}. 
For ~ iy\ ^ 1; ^(^) reads 

n{r) — >• -To 1;^ -T, as IJ — >• oo, 

e^'-l-e^ + e^' ' (3.23) 

n(0) = 2. 

From Eqs.(3.22) and (3.23), one can see that the angular momentum of the neutral vortex 



IS zero. 



We see from Eqs(3.9),(3.14),(3.16) and (3.20) that the charge and the angular mo- 



mentum Ly of our vortex have a relation : 



eLy = Qy, (3.24) 

in the large B limit. 



§4. Fluctuation corrections 

In this section we study the fluctuation corrections to the mean fleld solutions. We 
approximate the integration over U in Eq.(2.4) into a sum over the mean field solutions of 
Eq.(2.8) and the quadratic fluctuations around the mean fleld solutions. That is 

|^exp|-l(^)^^f/if/ + 0(if/^)l (4.1) 



n 



11 



where J7q"'^'s are solutions of Eq.(2.8) and is given by 

\ / n 



V )n~6U{y',x';t')5U{x,y;t) 

- V{y' - x')Dn{y', x', t'; x, y, t)V{x - y) = V{1 - D^V), 



= 5(x' - x)5{y' - y)5(t' - t)V{x - y) 
u^u^^^ (4.2) 



Dn{y',x',t';x,y,t) = l^'ip\y']t')'ip{x']t')-il)\x]t)'ip{y]t)\ 

Vt(^'; t')i;{x'; t')) U^ix; t)i^{y- 1)) (4.3) 



= (V't(y';OV'(y;i))^(„, {i^{x';t')^p\x;t)^^^^^ . 
The functional Gaussian integration in Eq.(4.1) is calculated as 



1 11 

exp[-- Tr ln(l - DnV)] = exp[- ^ - Tr(Z}„\/)^], 

(4.4) 



where 



1 If 

— TY{DnVy = — dti(fxi(fyi ■ ■ ■ dtj(fxjSyjV{xi - yi)Dn{yi,xi,ti] X2, 2/2, ^2) , , 
J (4.5) 

• • • ^i^j - yj)Dn{yj: Xj.tj] xi, yi, ti). 

Especially the j — 1 term is the first order of V and only includes the equal time propagator 
as 

- Tr{DnV) = -j dtd^xd'y (^ij\y- t)ij{y- 1))^^„^ V{x - y) (^(x; t)ij\x- 1)) ^^^^ 



^-^J d'xd^ut\x,x)V{x-y)ut\y,y). 



(4.6) 



Thus the j = 1 term is a direct interaction term. Graphical representations of Eqs.(4.5) and 
(4.6) are given in Fig.l. The energy eigenvalue appears in the partition function as e~^^" 
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for /3 — > oo. Thus 



£;tot =£;mean ^ ^ex ^ ^dir ^ ^^^2)^ 

M 



(4.7) 

Et dVy{U^''\x,x) - Uo)V{x - y){U^''\y,y) - Uo). 

is an exchange energy and E^^^ is a direct energy. In E^^^, we replace Uq"'\x,x) by 

(n) (n) 

Uq (x, x) — Uq, because the electron density Uq {x,x) becomes a uniform density Uq at 
infinity and V{r) does not couple the zero mode of the Fourier component. 

From a dimensional analysis, term behaves as B^^^^"^ and J = 1 term is a leading 
term of the perturbation expansion in the large B limit. 

§5. Numerical results 

We will solve the self-consistency condition Eq.(2.8) numerically for the charged and 
neutral vortex states of Eqs.(3.8),(3.15) and (3.21). At the start we expand F(p^) in Eq.(2.12) 
around — eB which is the energy of the lowest Landau level and approximate as 

= FieB) + F'{eB)ip'' - eB) = Fq - 
7r3/2 

Fo = e^C/o— — X 0.88934, (5.1) 

J_ ^ ^Uoi-f/^ X 0.48861. 
2m' 8k "V 

In this approximation, the eigenvalue equation in Eq.(2.13) becomes 

r(P + eA)2 (p + eA + n)2 ~ 

^ 2^ ^ W ^ °^ ~ PW)]a;-.a;'«Z,o(2;) = EiUi^o{x), 

uifii^) = vifl{r)e~'^^ ; Z = integer, 
where we add a constant Fq to the energy /io for a simplicity and put Ei — Ei q + Fq. For 
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the energy eigenvalue of the uniform density state, we define El''' as 

where M — mm! /{m + m'). We use the analytic form of Eqs.(3.12),(3.19) and (3.23) for 
ni{x) and solve Eq.(5.2) self-consistently with respect to p. We restrict the form of p{x,y) 
as 

p{x,y) = f{^/rvr^,0'^), 

(5.4) 

/(r,0) = p(r), 
and calculate p by a numerical self-consistent method 

(A) Charged vortex for riy = 1,2 

Using Eqs.(3.12) and (3.19), we obtain the numerical solutions of p. See Figs.2 and 3. 
The electron densities for the quasihole and quasiparticle vanish at the vortex core. Figures 
5 and 6 show the energy eigenvalues Ei. For large I, Ei becomes degenerate to E^^\ For 
I < —1, El belongs to higher energy levels and is larger than Ei for Z > by 0{eB/M). 
The energy differences from the uniform density state, AEn = E^ — E^^^ = E^ — E^^\ are 
listed in Table I for = 10 [T]. The gap energy A is a half of the pair creation energy of a 
quasihole and quasiparticle which are separated infinitely as 

A^/,^l{AE'^% + AE'_^l/,). (5.5) 

Experimentally A is measured as an activation energy by the temperature dependence of 
the diagonal resistivity 3S Pxx '-^ 

e ^ . We calculate A^^y3(the unit is kelvin) at i? = 
5, 10, 15, 20 [T]. See Fig. 8. The gap energy A2/3 is about two times as large as A 1/3. 

(B) Neutral vortex for Uy = 2 

Using Eq.(3.23) we obtain the numerical solution of p. See Fig. 4. Figure 7 shows the 
energy eigenvalue Ei. The energy differences from the uniform density state arc listed in 
Table.I for B — 10[T]. The gap energy Aneu for a neutral vortex is equal to AE^^. See 
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Fig. 8. It seems that our neutral vortex(charge = 0, 71^ — 2) is a bound state of two charged 
vortices (charge = — e/3 + e/3, nv = l + l). The gap energy Aneu is about a half of the pair 
creation energy 2A]^/3. Since the sum of energies for two charged vortices is larger than the 
energy of one neutral vortex (A£J^°|^3 + AE^j^^g > A£'*g^), it is possible to make a bound 
state. 



§6. Summary and discussion 

In this paper we applied a mean field theory to the FQHE for i/ = 1/3 and obtained the 
charged and neutral vortex solutions numerically. These solutions have fractional charges 
and fractional angular momenta and have a relation (3.24). We calculated the gap energies 
for these vortices. The smallest one among the gap energies is A for the pair of vortices 
which have charges ±e/3. The B dependence of A1/3 is written as 

Ai/3 = C~V^, (6.1) 

and C ~ 0.09. This result is consistent with other theoretical calculations (C ~ 0.1) based 
on the Laughlin's theory ^^'^^ and on the exact diagonalization of the finite system-^^^'^-*^). 
The experimental values for the activation energy are about half of these theoretical values 
at B — 10 ~ 20[T']. This discrepancy is explained by the effect of a thickness of 2D layer 
and effect of disorder Below B ~ 5[T] the experimental gap energy becomes very 
small and the discrepancy still remains. In Ref.5), the gap energy vanished at i? = 5.5[T] 
without the thickness effect and disorder effect. However we obtained negative results in 
the present paper by calculating fluctuation corrections. The gap energy did not vanish at 
B = 5 ~ 20[T]. 

In addition to charged vortices we obtained a neutral vortex which has a zero charge 
and a zero angular momentum. The B dependence of Aneu is written as Eq.(6.1) and 
C'neu ~ 0.11. The neutral excitation in the FQHE was studied by means of the single-mode 

15 



approximation and the exact diagonalization of the finite system These studies show- 
that the excitation energy has a magnetoroton minimum at finite wave vector (Croton ~ 0.15). 
We speculate that our neutral vortex may be a wave packet made out of the collective mode 
around the roton minimum. Recently the q — collective mode in the FQHE was observed 
by inelastic fight scattering^^) . In Ref.9) it is speculated that the q = excitation may be a 
two-roton bound state. We hope that the magnetoroton will be observed by inelastic light 
scattering or other methods. 
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Appendix 

We see in §5 that the energy eigenvalues become non-degenerate for small I. These cause 
a difficulty that many-electron state l^*) is not the eigenstate of Hjjo in Eq.(2.12). We can 
overcome this difficulty by using an arbitrariness in the limiting procedure of discretized 
time slices. In this appendix we mainly use the notation of Ref.6). 

The partition function can be written as 

Z = Tr(e-^^) = lim Tr[l - ip^Kip - -ip^^Vipip]^ , (A.l) 
where /? — Ne and we use the shorthand notations : 

^2 



ip^ili^V'ip'ip — / d^xd^y'^\x)'tp\y)V{x — y)ip{y)'4^{x). 



(A.2) 



Using an auxiliary field a, Z becomes 

N 



Z = lim / ^e~iEfc Tr[fT(l - eV^^^^ + eV^V^vV^ - — (TfcV^^V^VV^V^afc)], 
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where 

tp^akvip^ J (fx(fyii)^x)ak{x,y)v{x,y)'ilj{y), (A.4) 

ak'ip'^'ip^Vip'ipak = J d'^xd'^y(fx'(fy'ak{x',y')'ip^{x')'ip^{y)V{x - y)'ip{y')ip{x)ak{y,x). 

The suffix k specifies a discrete time slice. We can choose a function v{x,y) in Eq.(A.3) 
arbitrarily because the term Tjj^akVijj is linear in a and vanishes after integrating over a. For 
a general v we cannot exchange the order of the limit e ^ and integration f Da. Therefore 
we have to integrate over a for a finite e and take a limit e ^ in the end. With a fixed a 
the e^-term seems not to contribute to Eq.(A.3). However the integration with a Gaussian 
factor e~§^'^*'^'= makes a become of order e~^. Thus the e^-term becomes of order e by 
integration J Da and is not negligible. We change the integral variable as 

(^k{x,y) = a^^\x,y) +^k{x,y), (A.5) 

where a^^^ is a time-independent mean field solution which is determined in the following. 
The mean field (7^°^ is order e° and is order e~^, then is also order e~^. Hence, Z 



becomes 



k=l 



(A.6) 



= lim /^e-^-(-'°^'^). 
a^^^ is fixed by the stationary condition : 



In the limit e — > 0, cr*-^'* is given by 



= 0. (A.7) 



= v{y,x) l^ilj\y)'ilj{x)) ^^^^ , (A.8) 
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H^{o) = j (fxtp'^Ktp - J d^x(fytp\x) (i^\y)ip{x)'^ ^^^^ v{x,y)v{y,x)^p{x) (A.9) 
+ ^ / d'^xd^yv{x,y)v{y,x) (ip\y)ijj{x)) (ip\x)ip{y) 



We choose v{x, y) as 



M 

v{x, y)v{y, x) = V{x-y) + J2 ^^^W,o(l/)«/,o(^-) / (^^(?/)^(^))^(o) " ^"^"^O) 

1=0 



Using the ansatz form (2.10) for {^\y)%l){x)) , the energy eigenvalue equation for a single- 
particle becomes 

[K - F{{P+a)^)p{x,x) - AEiflU^^ui^oix) = El^uiflix), (A.ll) 
where we choosed AEi^ as 

AEi,o = Ei,o-eI'^. (A. 12) 

E^^Q is the eigenvalue for a uniform density state and is independent of I. Thus the energy 
eigenvalues become degenerate in / and the state |^) of Eq.(2.14) becomes eigenstate of H^m 
of Eq.(A.9). Since a difference between the eigenvalue equation in Eqs.(2.13) and (A.ll) is 
only a constant A.Eifi, the eigenfunction uifi{x) is common in Eqs.(2.13) and (A.ll). Prom 
Eqs.(A.lO) and (A.ll), the mean field energy and the exchange energy become 

M 

ET = \j d\d^yv{x,y)v{y,x) (i^\y)i^{x)) ^^^^ {i,\x)i,{y)) ^^^^ (A.13) 
= - d^xd^yV{x-y)U^^\x,y)U^''\y,x) + -Y,^AEi^^. 

1=0 
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The quadratic fluctuation around (T*^°) is calculated as 



^Skk'S{x' - x)5{y' -y)- e(l - 5uk')v{y' , x')Dn{y' , x' , t'^, x, y, tk)v{x, y) 
+ e8kk'S{y',x']x,y), 



S{y', x'; X, y) = U\y')^l^\x)V {x - ymx')^{y) 



(0) 



+ v{y',x')v{x,y) (ilj\y')ilj{x')J ^^^(^ij\x)i}{y) ^ 

(A.14) 

Prom Eq.(A.14), the direct energy included in the fluctuation correction is written as 

fi r fi (A.15) 

= 1 / cfxdMut\^,^)-Uo)V{x-y){U^''\y,y)-Uo) + ^Y.''l^^lfi- 

1=0 

Prom Eqs.(A.13) and (A.15), the total energy becomes 

Thus the results obtained in this paper are correct for the total energy E^°^. 
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FIGURE CAPTIONS 

1) (a) Feynman diagram for Eq.(4.5), j — 5. 
(b) Feynman diagram for Eq.(4.6), j — 1. 

2) P-n^/s for quasiholes oi Uy — 1,2 a,t B — 10[T]. 

3) p+n^/s for quasiparticles of riy = 1, 2 at = 10[T]. 

4) Pneu for a neutral vortex of ny = 2 at S = 10 [T]. 

5) El for quasiholes of = 1, 2 at S = 10[T]. o represents an eigenvalue for riy — 1 and 
• for Uy — 2. 

6) El for quasiparticles oi riy — 1,2 sit B — 10[T]. o represents an eigenvalue for Uy — 1 
and • for Uy — 2. 

7) £'^^0 for a neutral vortex oi Uy — 2 a,t B — 10[T]. 

8) The gap energy A^^^ for a charged vortex pair of ny = 1,2 and Aneu for a neutral 
vortex at S = 5, 10, 15, 20[T]. 
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